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We study the out of equilibrium dynamics of the infinite range quantum Heisenberg spin 
glass model coupled to a thermal relaxation bath. The SU(2) spin algebra is generalized to 
SU{N) and we analyze the large- A'^ limit. The model displays a dynamical phase transition 
between a paramagnetic and a glassy phase. In the latter, the system remains out of equilibrium 
and displays an aging phenomenon, which we characterize using both analytical and numerical 
methods. In the aging regime, the quantum fluctuation-dissipation relation is violated and 
replaced at very long time by its classical generalization, as in models involving simple spin 
, algebras studied previously. We also discuss the effect of a flnite coupling to the relaxation 

baths and their possible forms. This work completes and justifies previous studies on this 
model using a static approach. 
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The study of the non-equilibrium dynamics of classical glassy systems has been the subject of an intense research in 
the last decade. A lot of progress has been made using scaling arguments, phenomenological approaches and mean 
field theory. One of the major achievement is the theoretical explanation of the aging phenomena, which is one of the 
, most striking feature of glassy systems. The analysis of the out of equilibrium of (classical) mean field spin glasses 
has played a major role for several reasons. It has furnished a framework to understand, interpret and analyze the 
r \ experimental results and it has given important predictions on the violation and the generalization of the fluctuation 
' O ■ dissipation relation out of equilibrium |^ which has been experimentally tested recently j^] . 

Usually, many glassy systems can be analyzed within a classical approach since they are characterized by transition 
temperatures at which quantum mechanical effects are not relevant. Nevertheless, there are also interesting cases in 
which the critical temperature can be lowered to zero tuning a parameter which controls the strength of quantum 
fluctuations. This give rise to a quantum critical point at zero temperature 0|. Close to this point, the quantum 
fluctuations are very important and cannot be neglected. One example which has received much attention recently 
^ , is the insulating magnetic compound LiIl0a;Yi_3:F4 which is an experimental realization of an Ising spin glass in 
' a transverse field Other systems where glassy properties in the presence of quantum fluctuations have been 
] observed are mixed hydrogen bonded ferro-antiferro electric crystals Q , interacting electron systems Q , cuprates like 
lO ' La2-a;Sr3;Cu04 amorphous insulators 

' The theoretical study of quantum glassy systems has been performed following two different and complementary 
, routes. One dimensional models (like the Random Transverse Ising spin chain) has been extensively studied and it 
' has been shown that the Griffiths-McCoy singularities are very important close to the quantum critical point [ p^ . 
On the other hand, after the work of Bray and Moore pT| much attention has been focused on infinite dimensional 
(mean field) models |l|, |l|, |l|, ^ ITfll, 11. H H- In particular, recently, it has been shown 

p7[ that for the quantum spherical p-spin glass model the quantum fluctuations drive the transition toward a first 
order quantum phase transition at low temperature. The same phenomenon has been observed experimentally for the 
(~| insulating magnetic compound LiHoa;Yi_2;F4 In ]28t it has been argued that this phenomenon is to be expected 
O in a large class of systems. 

^ In contrast, the study of real time out-of-equilibrium dynamics of quantum glassy system is a recent subject 

and only very few results are available at the time of this writing. In a flrst pioneering paper, Cugliandolo and 
Lozano p5[ presented a detailed solution of a quantum version of the p— spin model. They showed how the out of 
. , equilibrium behavior of classical glassy systems is affected by quantum fluctuations. In particular they found that the 
■ low temperature glassy phase is characterized by the aging phenomenon. In this regime, the fluctuation dissipation 
relation is violated and it is generalized to a form that coincides with the (generalized) classical one. This could 
seem natural since at low frequency the quantum fluctuation relation coincides with its 7i ^ limit (for a bosonic 
system). Indeed it has been shown in ^ that for models with simple commutation relations (particles and rotors) 
the classical nature of the generalized fluctuation dissipation relation is due to the fact the dynamical equations are 
flxed point of the re-parameterization group of time transformations and the renormalized aging dynamics becomes 
classical at the flxed point. The quantum mechanics enters only as a renormalization of the coefficients of the 
dynamical equations. 

However what happens for models with a non trivial spin algebras, as the SU{N) model studied in ||2^, remained 
an open question. The study of the out of equilibrium dynamics of this type of quantum glassy systems is the main 
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aim of this paper. We will focus on the quantum Heisenberg Spin Glass where the SU{2) spin symmetry group is 
replaced by SU{N) and take the large A'^-limit. In this model, the spin are true quantum spins, i.e. with non trivial 
commutation relations, and this introduces in the problem Berry phases which play an important role Recently a 
detailed mean field solution using an equilibrium approach has been presented in [2^ . The model displays a second 
order phase transition at a temperature T^q between a paramagnetic phase and a spin glass phase, and it is solved by 
a one-step replica symmetry breaking scheme. Moreover, using a procedure called "the marginality condition" , the 
existence of a dynamical transition has been predicted at a temperature > Tf,q . First introduced and discussed in 
the quantum case in , this prescription was used in ||2|, ^ since it lead to the most acceptable solutions. Recently, 
the TAP approach has been fully generalized to quantum systems The relationship between TAP and replica 

approaches gives a further hint on why one has to choose the marginal solution in the replica method. In fact this 
solution is related to the marginally stable TAP states which have some flat directions around them in the (quantum) 
free energy landscape, on the contrary of all the others which are completely stable. Assuming that the quantum out 
of equilibrium dynamics is dominated by the presence of flat directions around the marginally stable TAP states, as 
it happens in the classical case, one finds a more natural justification of "the marginality condition" . However, only 
a complete dynamical analysis can fully justify this procedure. The analysis performed in this paper of the real time 
out of equilibrium dynamics, using the Schwinger-Keldysh close time formalism [ |3l| , shows indeed its correctness. 

This paper is organized as follows : in section ||, we present the model and the relaxation bath coupled to it. 
In section we present the dynamical large- iV equations for the retarded and Keldysh correlation functions and 
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we explain their derivation and how to deduce them from the simpler imaginary time equations. In section III 
present both an analytical and a numerical analysis of the dynamical equations. Numerical evidence for aging and for 
a generalized fluctuation-dissipation theorem in the aging regime are presented. Moreover the analysis of the aging 



regime justify the "marginality condition" used in previous works pl| , [22| . Finally, in section IV, we briefly discuss 
the effect of a finite coupling to the relaxation bath. 

I. THE HEISENBERG SPIN GLASS AND THE RELAXATION BATH 

The model considered in this paper is a quantum Heisenberg spin-glass on a completely connected lattice of A/" sites 
with quenched disordered couplings Jij / V MN which are independent random Gaussian variable of zero mean and a 
variance Jfj/{AfN) = j'jj/iMN). Each spin is linearly coupled to a thermal bath. Moreover we generalize the SU{2) 
spin symmetry group to SU{N) and we take the large A^- limit. This generalization allows us to obtain a tractable 
model which has still highly non trivial quantum effects and reproduces qualitatively well the known results for the 
SU{2) model, as far as it has been possible to compare the N = 2 and the N = oo cases [p^ [2l| , The Hamiltonian 
reads: 

^ " ^ "^'^ ' ^ ^ ^' ' ^"^ ^ 

where the scaling of the spin-spin couplings and the antiferromagnetic spin-bath coupling has been chosen in such a 
way to obtain a sensible large Af, N limit, i.e. Jfj, Jb oc 0{1). The first term is the Quantum Heisenberg spin glass 
Hamiltonian, the third and the second terms represent respectively the thermal bath of spins Si and its coupling to 
the spins Si via the coupling constant Jb- Let us now discuss them separately. 

Among the possible representation of the SU (N) spin, two versions have been studied ^ : the bosonic model, 
in which the spin operator S is represented using constrained Schwinger bosons b by 

Sat3 = bl^bfj - SSa[3 (2a) 

N 

bib^ = SN (2b) 

a=l 

and the fermionic model, in which the spin operator S is represented similarly using Abrikosov fermions / by Sa[3 = 
ftfl3~ ^^a0j "^ith the constraint J2a fafa = SN (0 < S < 1). See Q for an introduction to this two representations. 
The two models are technically very similar but there is an important physical difference between them : in the 
fermionic model, quantum fluctuations are so strong in the large- A^ limit that the spin glass ordering is destroyed [3^ 
(the critical temperature vanishes when A'^ diverges p^ ), whereas in the bosonic model, a spin glass phase exists at 
low temperature [pT|, p3] . In the following we will focus mainly on the latter one and we briefly discuss some results 



for the former one at the end of Section HI. In the model we study, the size S of the spin is a fixed, tunable parameter 



which controls the strength of quantum fluctuations (since AT = oo it is moreover continuous). Let us emphasize that 



3 



the N ^ ao limit is not the classical limit : as shown in ^ , the model is classical for large S (but not a very low 
temperature), while it is "more quantum" (i.e. the quantum fluctuations are more important) at low S and displays 
a quantum critical point at S* = where the spin glass temperature vanishes. 

Moreover, it is important to remark that the real physical object is the spin S, not to the boson b or the fermion / 
which should be considered here more as mathematical tools. Technically, this leads to a simple U{1) gauge invariance 
of the bosonic or fcrmionic theory (we can always multiply the b or the / by a phase) whose consequences will be 
explained later. 

Let us now discuss the role of the relaxation bath terms in (^. Its presence is necessary to allow the energy 
dissipation. It guarantees the relaxation toward equilibrium above the dynamical transition Td and it is required to 
obtain an aging regime below Tj^. In this paper, we are mostly interested in the Jb limit, which must always 
be taken after the long time limit : for example, the dependance of the equilibrium state in Jg is expected to be 
smooth in this limit, although the transient time towards the equilibrium diverges. The bath we have considered in 
(|l|) is supposed, as usual, to be very big and always in equilibrium at a finite temperature T = 1//3. For further 
simplification we take independent baths from site to site (labeled by i), and the spins Si carry an additional degree 
of freedom a, with 1 < a < N"f (where 7 is a constant), which ensures that the bath is much bigger than the spins 
it is coupled to. Moreover we will make the assumption of factorized initial condition pst . By this we mean that the 
initial density matrix is a product of an equilibrium density matrix for the bath and an initial density matrix for the 
system. Since 6 (or /) is not the physical object, the bath should respect the U{1) invariance, i.e. it must couple 
to two 6's and we do not consider baths coupling linearly to b in this paper. Of course many different choices are 
possible. We first consider a "generic" bath Hbath{s) of interacting spins s and expand the Keldysh effective action for 
S at second order in the coupling constant Jb (the first order vanishes). In the dynamics, the bath will then appear 
only through its susceptibility xo (see section |l|). This approach is appropriate when we only consider the bath as a 
device to provide thermalization. One could wonder how correct is to study a low temperature glassy phase using a 
perturbative treatment of the coupling to the environmental heat bath. In Section IV we will show that this is not a 
limitation and we will briefly discuss the simplest type of spin bath which will turns out to be a Kondo bath. 



II. THE DYNAMICAL EQUATIONS 

In order to study the real time dynamics, we use the Keldysh method : the time evolution operator is written as 
a path integral over 2 times i+ and t_ , running from to infinity forward and backward respectively. In the classical 
limit, this method reduces to the Martin-Siggia-Rose-DeDominicis-Janssen formalism (See Appendix C of |^). We 
take an infinite temperature initial condition at t = and do a instantaneous quench to temperature T. This means 
that the initial density matrix for the system is simply the identity operator. As a consequence the initial density 
matrix, which is a product of the equilibrium density matrix of the bath and the infinite temperature density matrix 
of the system, does not depends on the disordered couplings. Therefore, as in the classical case, there is no need to 
introduce replicas to compute disorder-averaged quantities. 

The quantities that we want to compute are the averaged response and correlation of the spin, which are defined 
as (the spins are on the same site) : 



Rs{t,t') 



^9{t - t')(s{t) ■ S{t') - S{t') ■ S{t)) = ^Y.(^it- t')([S^,{t), Sf^^t')]) 
CsM ^ :^{s{t) ■ S{t') + Sit') ■ Sit)) = ^ 5] {Sc.f3it)Sf,^it') + S^f3it')Si3o.it) 



(3a) 
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where the bar denotes the average over disorder and the brackets denote the "Keldysh average" i.e. the Hamiltonian 
evolution of the quantity starting from the initial condition at i = 0. Moreover, in this paper, we take h = 1. 

In the following, we first derive the Keldysh action, we average over disorder and we take the large-iV limit; then 
we recall the so-called "Larkin-Ovchinnikov representation" and we express the dynamical large- equations in their 
final form using the retarded and Keldysh functions of the bosons R and K from which one can obtain Rs and C5. 
It is not possible to obtain tractable equations for the physical quantities Rs and Cs directly, contrary to rotors ]3C| ] 
or spins models, and that makes the problem more complicated |2^ . 

We start using a Keldysh action defined on the double contour : 



Sr, S r 



dt 



E 



snt)-snt) 



Jb 



J2snt)-^jt) + H^M^it)) 



(4) 
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In this expression, a = ± denotes the upper/lower contour, Sg are the Berry phase on the upper and lower contours 
p6| . After the average over the disorder we take a saddle point over the number of sites Af. Hence, we get a 
self-consistent problem (some scalar products have been explicitly written with SU (N) indices for clarity) : 

E <^l°" dt^(^S^it)-^Jt)+H^,,^{^{t))^ (5) 

where i-)^^ means the average over the single site action (^. 

At this stage, it is useful to define the spin correlation functions in the so called "± representation" : 



x°.(M')^^(r.-«w..-*(0)B.th (7) 

where a, 6 = ± are the contour indices, and T is the time ordering on the double contour, {■)-qi^^yi denotes the average 
with respect to the bath. Now, using explicitly the SU{N) invariance of the theory, integrating out the bath, and 
expanding to second order in Jg we get an action only for the spins (summation over a and /? is implicit) : 

^E«^ ir dtdt'Slp{t)Sl^{t')[jlxlb{t:t^ (8) 



a.b 



Until now, the derivation is correct for any value of N and in particular for N = 2. The great technical advantage of 
the large- TV limit becomes manifest if one considers in detail the self-consistent single site problem. In fact, because of 
the presence of the Berry phase, the single site measure defined by the action (^ is far from being simple. Indeed the 
single-site functional integral cannot be performed and as a consequence it is not possible to obtain a closed equation 
for the spin-spin correlation function. The large- iV limit simplifies the single-site measure and gives a set of closed 
equation on the two-point functions. Using the Schwinger bosons, the Berry phase contribution to (|^) reads: 

SB = -Y.a dtb-Jdtb-^ (9) 

while the other part of the action can be obtained simply replacing with its expression in terms of bosons. For 
a finite TV the problems remains still very complicated since one has to integrate only on bosonic fields respecting 
the constraint ([2b|). Whereas in the large- limit, which we shall study in the following, the sum bl^ba/N does 
not fluctuate and this greatly simplifies the analysis. In particular we can obtain the saddle point equations on the 
bosonic Green functions, which in the ± representation are defined as : 

Gab{t,t') = -l(Tba{t)bl{t'))^ (10) 

The computation can be done explicitly, using the same decoupling as in the imaginary time equilibrium computation, 
as explained in [2^ . Here we present a faster derivation, noting that the same diagrams for the self energy derived 
in the imaginary time computation appears in the Keldysh formalism : a first diagram, corresponding to the spin 
glass interaction itself, and a second one, corresponding to the coupling to the relaxation bath (See Figure |l]). Thus 
using Feynman rules in the ± representation, we find immediately (the factors can be checked against Matsubara 
computations : see appendix A) : 

^abit, t') = - J^abGl,{t, t')Gba{t', t) + Jlabxlbit, t')Gab{t, t') (11) 

We see that the bath only enters though its susceptibility x°. For simplification, in the following, we take a 
specific form for the susceptibility of the bath Xabi^'^') ~ ^ahC^i ^')^ba(^'' where G° is the Green function of free 
fermions with a Lorentzian density of states at half filling (See section IV for a discussion on the relaxation bath and 
a justification of this formula). 



5 



^ab{t, t') = — J' 



Gab{t, t') 




+ 



Gab{t,t') 




a, t 



b,t' 



FIG. 1: Feynman diagrams for the self energy : tie solid line is the bosonic Green function, the dashed lines represents the 
x" of the bath (we take a product of two fermionic functions, see text). 



To simplify the analysis of the dynamical equations, it is useful to write them in a different way, using the so-called 
"Larkin-Ovchinnikov representation" (LO) of the equations, in which the (matrix) Green function is given by 



G = 



1 /I 1 



2 VI -1 



G+- 
G— 



1 1 

-1 1 



R K 
A 



(12) 



where the retarded (response), advanced and Keldysh (correlation) two-times Green functions are defined by : 

R{t, t') = -i0{t -t'){ [b{t), {t')] ) = G++ {t, t') - G+- {t, t') (13a) 
A{t, t') = +i9{t' - t) { [b{t), b\t')] ) = G++(t t') - G_+(t, t') (13b) 
K{t,t') = -t{{b{t),b^{t')})^G++{t,t') + G.^{t,t') (13c) 

Note that our convention for R differs from the one used for Rs- Moreover, we will use the relations : 



A{t,t') 



R{t',ty 
~K{t\ty 



(14) 
(15) 



to eliminate A and restrict ourselves to i > t' . The (LO) representation is simpler because it uses the relation 

G++ +G__ = G y + G^ to reduce the number of functions (just R and after that A is eliminated using (14) and 

because it makes the causality of the equations explicit. In the (LO) representation the Keldysh indices structure of 
the vertices are particularly simple : the 2-leg vertex is Sab , which leads to a simple Dyson equation G^^ — Gq "'^ — S 
(where the inverse are just matricial inverses), and the four-leg vertex used in ( pl| ) is {1 ® (Tx + cr^ ® l)/2, where 1 is 
the 2x2 identity matrix and is the usual Pauli matrix. Remarkably, this vertex factor is fully symmetric in the 
Keldysh space. Hence, a quick way to switch to the (LO) representation is to recompute the Feynman diagrams (See 
Figure (|l|)) with the (LO) Feynman rules (although a direct computation using just the definitions is possible but 
more tedious). After these manipulations, we finally obtain the main equations of our paper (for t > t') : 



The Dyson equations ((16c) is rewritten to involve only functions for t > t') 



idtR{t,t') = S{t -t') + J duT.B.{t,u)R{u,t') 

idtK{t,t')^ [ du j:ji{t,u)K{u,t') + [ du T,K{t,u)R* (t' ,u) 
Jq Jo 

du^R{t,u)K{u,t') + [ du (i:K{t,u)R*{t' ,u) -i:R{t,u)K*{t' ,u) 



(16a) 
(16b) 
(16c) 



The boundary conditions, which derive respectively from (2b) and the commutation relations of the boson 



K{t, t) ^ -i{2S+l) 



lim R(t,t') = -i 
*^(*')+ 



(16d) 
(16e) 
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• The self-energy in LO representation (these formulas are local in times, so the argument (t, t') has been omitted 
for clarity) : 



i:^^''' = -^\^[\^RY-\KY'jR+[R*K-K*R)K] (16g) 
Y}^^''^ = -^(^{\r\^-\k\^^K+[r*K-K*r]r] (16h) 

j^Bath ^3.((\r^\^^\k^\^\r+Ir*^k^^KMK] (16i) 



I 2 I I 2 



(16j) 



where i?o and are defined in the same way that R and K starting from f,. This expression emphasizes the very 
similar structure of the two terms in the self-energy : more generally, the bath term reads E^'**'' cx ~x!kR + X^rK and 
j^Bath ^ y^j^R -\- x'rKj where x'r and x^k the retarded and the Keldysh part of respectively. 

Finally, we now derive the expression for the response and correlation functions for spins defined in (pa|). In the 
iV — !■ CX) limit, Cs and Rs can be easily computed from R and K using the relations : 

a/3 

^lin{K{t,t')R*{t,t')) (17a) 

a/3 

= ^{\K{t,t'r-\Rit,t')f-\R{t',t)f) (17b) 



To derive (17a), we used Eqs. (Ea,pq), the SU{N) invariance, the limit N oo (to drop subdominant terms), and 



the relations between Ga.b and R,K,A that invert ( 13a 



The equations (y^ describe the dynamics of the Quantum Heisenberg spin glass in the large- iV limit, coupled to 
the bath. In the following sections, we present an analysis of these equations both in the paramagnetic regime and in 
the aging regime, together with some results extracted from a numerical solution of this systems. Let us first make a 
few preliminary remarks : 

i) Contrary to the quantum p— spin problem studied in [ p5[ , there is no need here for a Lagrange parameter 
associated to the constraint. This is due to the fact that the constraint (|2b|) commutes with t he Ham iltonian 
and hence is conserved in the time evolution. Indeed one can check that ( |16g 16h,16i 16] ) and (16a, 16c) imply: 



^^^cxIm|/^(t,OI,,=, = (18) 

Thus if the boundary condition is verified for t = t' = 0, it will propagate at all later times. 

Formally, one can introduce such a parameter A in the equation, by replacing idt by idt + X^ but it can be removed 
using a C/(l) gauge transformation [K{t,t'), R{t,t')] [K{t,t'), R{t,t')] = [e~''^^*-^'^ K{t,t'),e-'^^*-^'^ K{t,t')] 
: '\i K,R are a solution of the equations with \\K,R\ are a solution with A = 0. This symmetry comes from 
the fact that we represented the spin (the physical object) with the bosons (a mathematical tool) and that the 
bath couples to the spin and not to the boson, and thus can not break the symmetry; this has an important 
consequence (see ii). 

n) In equilibrium, the spin response and correlation functions are related by the quantum fluctuation dissipation 
relation QFDR: 

Rs{t) = tOir) ^ exp(-ic.r) tanh Cs{t) (19) 
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It is important to notice that instead the boson response and correlation functions R and K are not in principle 
related by this QFDR. Technically, this is due to the U{1) invariance explained above : if [K{t,t'), R{t,t')] 
satisfies QFDR [K, R] will not in general. Imposing QFDR for the boson demands that A takes a precise 
value Aq. However, this is not a problem since the only physical objects are the spin response and correlation 
functions. When using the Matsubara imaginary time formalism, one does not face this difficulty, since one 
automatically requires the QFDR to be satisfied, because of the /3— periodicity of the imaginary time boson Green 
function. Thus, when we will analytically continue our equations in imaginary time to compare to Matsubara 
computations (See Appendix ^ ), we will have to reintroduce Aq. 

in) The presence of the thermal bath is clearly required : for Jb = 0, the solution of the equations has the property 
K{t,t' ) = (25*+ l)R{t , t') for all t,t' (this can be check order by order in the coupling constants Jh and Jb, 
using (16a,16c,16g, r6q ,16i,16j)) and it is clearly incorrect (for example, it can not satisfy the high-temperature 
limit of the fluctuation-dissipation relation). 



Jjj {S is 



iv) We can immediately generalize these equations in the fermionic case by changing S 

the size of the "fermionic" spin, and the sign change in front of the coupling constant comes from the fact that 
there is now a fermion loop in the diagram). We will see in Section [II that this simple change leads to the 
disappearance of the aging phenomena, as expected pij . 



III. SOLUTIONS OF THE DYNAMICAL EQUATIONS 



In this section, we present the solution of the dynamical equations ( [iq ) using both numerical and analytical results. 
Indeed, these integro-differential equations are causal, so one can construct the solution step by step in time. This 
property is very general (See p5| for another example) and it is the basis for the numerical algorithms, although in this 
problem some new technical refinements are needed in order to compute an accurate solution at a reasonable cost in 
computational time (see the appendix ^ for a detailed discussion) . The numerical solution shows that the model has 
a dynamical phase transition at a temperature Td{S, Jb) between a paramagnetic phase (T > Td) and a glassy phase 
(T < Td), as expected on general grounds and predicted in [ pT| p2| . At high temperature, the system equilibrates 
inside the paramagnetic state: after a transient time all the two-time quantities become time-translation invariant 
(TTI) and the quantum fluctuation dissipation relation QFDR holds. Instead at low temperature the system never 
equilibrates on finite timescales |]37|| and one can identifies two different time sectors on which the two-time functions 
evolve: when t' is large but the difference t — t' is of the order of one and very small compared to t' the system 
seems to be equilibrated (the QFDR is approximatively verified and all the one time quantities, as the energy or 
the Edwards- Anderson parameter, have almost converged to their asymptotic values). However, on larger timescales, 
when t — t' is of the same order of t' an extremely slow dynamics sets in. In this regime the QFDR is violated 
and the aging phenomenon appears ||l], ^, We remark that if one takes the long time limit and then sets the 
coupling to the bath to zero then the dynamical solution gets back to the equilibrium for T > Td only. For T < Td 
the system never reaches a stationary solution. However, the pseudo-equilibrium solution reached in the time sector 
t — t' ^ 0{1) << t' can be also obtained by a pure static computation using the marginality prescription p^ . In 
the following, we will take Jb > {Jb = 1 for numerical computation). In section IV, we will discuss what happens 
changing the value of Jb- 



A. Equilibration into the paramagnetic state 



At high temperature, the numerical solution shows, as expected, that the system equilibrates into the paramagnetic 
state after a transient time t^q'- for t^t' » t^q the response and correlation becomes a function of r = t — t' on ly and 



they are related by the QFDR (see Fig. g). This is indeed what we obtain from the analysis of Eqs. ( 16a 16c ) in the 
limit t,t' —f 00 with t = t — t' fixed. In this limit the equation on K and R can be easily written in Fourier space (we 
reintroduce the A term, in agreement with the discussion at the end of section |ll|) : 

R^^{uj) = uj + \-Y.ii{uj) (20a) 
K{uj) = I]K(w)|i?(w)p (20b) 

As in [psf , it is possible to show order by order in perturbation theory that these equations admit a solution such 
that the spin correlation and response functions satisfy QFDR. In Figure we plot the spin correlation function 
Cs{tw + T,t.it,) and the response function Rs{tw + T,tj^) as a function of r for different for Js = 1, Jh = 1 
and T ~ 10. These figures represents the typical behavior of Cs and Rs in the paramagnetic phase: after a short 
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FIG. 2: Spin response as a function of t for = 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 10, 15, 20, 25 (from bottom to top (tmax ~ WO)), 
Jg — Ju = 1, S = 1 and T = 10 (inside the paramagnetic phase). After a transient time the response converges to a stationary 
equilibrimn regime. Inset: spin correlation as a function of r for the same parameters. 



transient time the functions become TTI (they do not depend on tu, anymore) and they decay quickly as a function 
of r. Moreover, in Figure ^, we plot the retarded function computed directly from the numerical solution and from 
the correlation function using the QDFR. The excellent agreement shows that QFDR is satisfied and the system has 
relaxed to equilibrium. 

This asymptotic solution represent the equilibrium dynamics inside the paramagnetic state and it exist only above T^. 
At T = Td the equilibration time diverges and remains infinite in all the low temperature phase (T < Td), as clearly 
indicated by the numerical solution (see Fig. ^). The study of this regime is the subject of the next subsections. 

In principle, one would like to take a small Jb, so that it allows the system to relax but does not change the value 
of the paramagnetic state. However, the relaxation time diverges when Jb goes to 0, even in the paramagnetic state, 
and this prevent the numerical program to converge towards the solution in a reasonable amount of time. We found 
that Jb = 1 is a good compromise. Indeed the imaginary time computation shows that the results at Js = are 
close to Js = 1. This relatively big value reflects the fact our bath, coupling to the spin degrees of freedom, is not 
very efficient. A more precise discussion will be given in Section [V. 



B. General properties of the glassy dynamics 

The numerical results and the analytical analysis of the dynamical equations indicates that the system remains 
always out of equilibrium at low temperature [T < Td). In the following we present the Ansatz which gives the 
asymptotic solution in the glassy regime and we compare it to the numerical results obtained integrating the dynamical 
equation numerically. This Ansatz is a slight generalization of the one introduced by Cugiiandolo and Lozano 
for quantum glassy systems, which is itself a generalization of the one discovered by Cugiiandolo and Kurchan for 
classical glassy systems 
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FIG. 3: Fourier transform of the spin response function (dots) compared to its expression (continuous line) computed by the 
QFDR from the correlation for Jb = Jh = 1, S = 1, tw > 10 and T = 10 (inside the paramagnetic phase). The excellent 
agreement shows that the system is fully equilibrated. 



The weak- erg odicity breaking and the weak long-term memory Ansatz 



In the long time limit (i,t' >> 1) we make the following Ansatz for the behavior of the bosonic correlation and 
response function pa]: 



K(t,t') 
R{t,t') 



Ksrit - t') + Kag 



(m 

{hit) 



-i\(t-t') 



lim KsT{t) = 

Kst{Q) = -i{2S+l) + ig 
Kag{0) = 



(21a) 

(21b) 

(21c) 

(21d) 
(21e) 
(21f) 



where h{t) is an increasing function of t, h'(t) is the first derivative of h{t), g is a real number |Q and 5^/4 is equal 
to the Edwards-Anderson parameter qea- Note the physics hidden in this Ansatz: in the time regime in which t,t' 
are large but their difference remains finite (called TTl-regime in the following) the system seems to have reached 
a stationary state, however on a timescale diverging with t,t' (t,t' are large but the ratio h(t')/h(t) remains finite) 
there is a secondary evolution called aging Q. This unveils the interpretation of h{t) as a self-generated effective time 
scale. 

Moreover we notice that if this scenario is realized for the bosonic correlation and response functions then it will be 
also realized for the spin correlation and response functions and for the self-energies. The presence of the oscillating 
exponential is a slight generalization with respect to and it could be gauged away, as discussed before. Moreover, 
when one compute the spin correlation and response functions these exponentials cancel. 

The second key ingredient that makes the asymptotic problem tractable is the assumption of the weak long-term 
memory property [|^, ^ that allows one to decouple the transitory regime from the asymptotic one. In fact the 
dynamical equations contains explicitly memory terms which couple all the timescales. So, how can one analyze 
the asymptotic regime without solving the complete problem? Within the weak long-term memory scenario the 
(linear) response to a finite time perturbation vanishes in the long time limit (t, t' — > cxd), whereas the response to a 
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perturbation which acts on infinite timescales (i.e. diverging as t, t') is finite. More precisely: 

hm / R{t,u)f{u)du^O but hm / R{t,u)f{u)du ^ , (22) 

where f{t) is a generic function. Therefore the dynamics on "infinite timescales" {t, t' —>■ oo) decouples from the 
transitory regime. 

Finally, we note that more general Ansatze with a set of different diverging timescales have been used in the context 
of classical |oj and quantum Q glassy systems. However, this type of solutions are physically and technically 
related to a full replica symmetry breaking solution in the thermodynamical analysis. For systems characterized by 
a one step replica symmetry breaking solution in the thermodynamics, as the model we are focusing on pl| , one 
generally expects only one diverging timescale. 



2 Generalized QFDR 



An outstanding physical property of the asymptotic dynamical solutions, discovered in the classical case by Cuglian- 
dolo and Kurchan ||^ and in the quantum case by Cugliandolo and Lozano |25|, is that in the aging time-sector {t, t' 
and t — t' are large but the ratio h{t')/h{t) stays finite) the standard fluctuation dissipation relation is violated but 
there exists a generalized fluctuation dissipation relation (GFDR) between the correlation and the response functions 
which has the usual functional form of the FDR (generalized to non TTI functions) and in which the temperature T 



is replaced with an effective temperature Tcff |25 . Two remarks are in order concerning Toff- First, a physical one: 
the effective temperature has a real physical meaning of temperature p| since is what a thermometer, whose reaction 
time equals the timescales on which the aging evolution takes place, would measure. Second, a technical one. The 
effective temperature T^g is related to the breaking point x arising in the replica symmetry breaking solution of the 
thermodynamics, i.e. Tcs = T/x. A general argument to show why one expects this to be true for a very large class 
of classical systems, included finite dimensional systems, has been presented in p2| . 

It is important to note, as pointed out in |Q that in the quantum case the GFDR is expected to become classical. The 
argument is the following: if Tog is finite then the Fourier integral relating the correlation and the response is dom- 
inated by w oc 0. Hence, one can develops the hyperbolic tangent recovering back a classical generalized fluctuation 
dissipation relation, which in our case (for spins) reads: 

Rs{t,t')^^dt'Cs{t,t') (23) 

J off 

This becomes a relation between the aging functions: 

Rs''{^^) = (24) 

where /i = h{t')/h{t). Moreover this has been argued to be generically true for models with simple commutation 
relations (particles and rotors) in ]30[| since the dynamical equations are fixed point of the re-parameterization group 
of time transformations and the renormalized aging dynamics becomes classical at the fixed point. The quantum 
mechanics enters only as a renormalization of the coefficients of the dynamical equations. We will show that this is 
also the case for our system which is characterized by non trivial commutation relations between the spins (contrary 
to the case of rotors or particles). However, the behavior next to the quantum critical point is still unclear (See section 

0)- 



C. Numerical results 

The numerical procedure used is described in Appendix It turns out that the problem is more difficult to solve 
than the classical ones or the quantum p— spins model, because the dynamical equations are for the auxiliary boson b 
and not directly for the physical spin S. As clearly indicated in the Ansatz, even in the aging time-regime where the 
spin correlation and response function evolve very slowly, the bosonic functions oscillate wildly. The numerical solution 
thus demands a more sophisticated algorithm, inspired from well-known methods to solve one variable differential 
equations. 

The numerical results support and validate completely the out of equilibrium scenario encoded in the Ansatz (IIIB 1). 
In the following we present the results for Jh = 1, Jb = 1 and a temperature well inside in the glassy phase T = 0.1. 
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FIG. 4: Spin correlation as a function of r for t^ = 15, 20, . . . , 95 (from bottom to top (tmax = WO)), Jb ~ Jh = S — 1 and 
T = 0.1 (inside the glassy phase). The height of the dotted straight line equals the Edwards-Anderson parameter computed 
within the static formalism and coincides well with the plateau value. The aging behavior is explicit. Inset: zoom of the spin 
correlation as a function of r on the time interval corresponding to the stationary regime for t^ ~ 5, 10, 25, 37.5, 50, 75 from 
bottom to top. The curves show a clear convergence toward a TTI stationary regime. 

In Fig. ^, we plot respectively the spin correlation function Cs{tw +T,t^) and the spin integrated response function 
xi^w +T,t^) = //""''^ dsRs{tw + s) as a function of r for different values of t^. We remark that a pseudo stationary 
regime sets in for r < 1 with a plateau, whose height is the Edwards- Anderson parameter. Its value (gl;*^ — 0.88), 
computed from the static analysis by the marginality prescription, is represented with a dashed line on Figure ^ : this 
shows a very good agreement between the two methods. 

Moreover the fact that correlation and response are related by the QFDR in this time sector, see Fig. shows 
nicely that this is indeed a pseudo-equilibrium regime. 

On a longer timescale t (X tw the aging behavior sets in as clearly indicated in the two figures. The integrated 
response and the correlation evolve more and more slowly increasing tw . Note that the behavior of x shows explicitly 
the weak long term memory scenario: even if the response vanish in the aging time sector the integrated response does 
not. Therefore the magnetization response to a constant and small magnetic field switched on at tw depends always 
explicitly on and becomes slower increasing . Moreover we note that the numerical results strongly suggest that 
h{t) ~ t. Indeed the different curves Cs{tw + T,tw) collapse very well on a single curve {Kag{iJ')) when plotted as a 
function of (t^, -t- T)/tw, see Fig. ^ We have also verified the correctness of the GFDR hypothesis. In Fig. ||, we show 
a parametric plot of x as a function of Cs for different values of t^- For classical systems [0 the limiting curve (for 

— > oo) is very useful to characterize the aging behavior. For systems with only one diverging timescale, because of 
the form of the FDR and the GFDR, one finds a two straight line plot where x c< — yC's for qsA < Cs < C{t, t) and 
X oc — Cs for < Cs < Qea- In the quantum case the situation is more involved since in the stationary regime 
Rs and Cs are related by the QFDR, therefore one does not expect that for qea < Cs < K(t,t) the plot should 
be very useful. However, as discussed in [^, in the aging time sector where < Cs < Qea the GFDR becomes 
classical and one should recover a straight line for x whose slope equals the inverse of the effective temperature. This 
is indeed what we find in Fig. |^, where we compare the behavior of x at low Cs with a straight line whose slope is 
l/7cff = x/T. X is the breakpoint in the one-step replica symmetry breaking solution obtained (for the same value 
of the parameters) generalizing the analysis performed in |2l[] to taking into account the presence of the bath (see 
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FIG. 5: Spin response as a function of r for t„ = 5, 10, 25, 37.5, 50, 75 (from top to bottom (tmax = 100 )), Jb ~ Jh = S — 1 and 
T = 0.1 (inside the glassy phase). Inset: zoom of the integrated response as a function of t on the time interval corresponding 
to the stationary regime for t^ = 5, 10, 25, 37.5, 50, 75 from bottom to top. The aging behavior and the weak long term memory 
scenario are explicit. 



Appendix ^). 

Finally, we can obtain the equations for the fermionic model by the simple change S — > ~S,J^ 
as discussed previously. Numerically solving these equations we have found no glassy behavior as predicted in [ p4[ . 
Indeed we show in Fig. ^ the spin correlation function in the fermionic case, which does not show any aging behavior 
at low temperature, whereas the same calculation for bosons (with the same parameters) clearly does. 



D. Analysis of the stationary regime 



We foc us now o n the time sector in which the difference between t and t' stays finite and i, f are very large. Hence, 
in Eqs. (21a, 21b) the aging part does not evolve and is zero for the response and equals —ig for the correlation. 
Plugging the Ansatz ( 2l421b[ ) into the dynamical equations we get: 



{id + X)RsTit) 
{id + X)KsT{t) 
A 



dt'j:%^{t-t')RsT{t') + s{t) 



/ + 00 f + oo 

dt'Y.j^(t^t')KsT{t')+ / dt'Y.i'^(t + t')R*sT{t') + A 



lim 

t — ^oo 



du 



h{t) " \h{t) J ^"^yhiu); t^oo 



dui:jP 



(h{u)\ h'{u) 



\ h{t) J h{t) 



-R 



AG 



h{u) \ 
h{t)) 



-ig dt'T.fi^{t') + j:^ dt'R^'^{t')+i\g 
Jo Jo 



(25) 
(26) 

) 

(27) 



where we have used for the self-energy the same notation introduced in (21a, 21l| ), and = E;^'^(l). Since we 
have defined Ksrit) in such a way that it vanishes in the long time limit, A has to be equal to zero. Note that this 
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FIG. 6: Parametric plot of the spin integrated response as a function of Cs for tu, — 40,45,55,60,65,70 (tmax = 100 j, 
= Ju = S = 1 and T = 0.1 (inside the glassy phase). The collapse predicted by the dynamical Ansatz is good. The vertical 
dotted straight line indicates the values of the Edwards-Anderson computed within the static formalism. The dashed straight 
line has a slope —l/Tcg where l/TcS — x/T and x has been computed within the static formalism. The curves clearly show 
that the generalization of the Huctuation dissipation relation holds in the aging regime. 

overall equation couples the stationary and the aging regime. As in the paramagnetic case, one can show (order by 
order in perturbation theory) that Rst and Kst satisfy the QFDR (and therefore Rg^,C§^ too). It could seem 
that the procedure to fix A is different in the two cases. In fact in the paramagnetic case one chooses A in such a way 
that the QFDR is verified for bosonic functions, instead now A is such that the correlation and response functions 
do not oscillate in the large time limit when t — t' and i, t' are very large. But the two values of A are the same an 
asymptotically oscillating function cannot satisfy the QFDR relation. 

As a conclusion the stationary equations can be fully interpreted as equilibrium dynamical equations. Indeed it 
has been shown in the classical ||^ and recently in the quantum case that this type of equations represents 
the pseudo-equilibrium relaxation inside the marginally stable TAP states (local minima of the free energy landscape 
whose Hessian is characterized by a vanishing fraction of zero modes). Imposing the marginality condition in the static 
computation is equivalent to consider a Boltzmann measure restricted to the marginally stable TAP states. It is for 
this reason that one can get information about the out of equilibrium dynamics by a purely equilibrium computation. 
However, it is important to understand that the equations ( |2^ , |2^ , ^ ) do not really represent an equilibrium relaxation. 
Because the marginally stable TAP states have a vanishing fraction of zero modes, the system find always a way to 
"escape" to these states, even if more and more slowly and this gives rise to the aging behavior. Hence, the physical 
mechanism inducing the slow dynamics is not an activated jump dynamics across some energy barriers but it is an 
entropic effect. The slow dynamics of the system is due to the fact that the longer is the time the smaller is the number 
of directions along which the system can escape. Finally, the fact that the marginally stable TAP states dominate 
the off-equilibrium dynamics whereas they are not relevant for equilibrium properties helps to understand why the 
dynamical transition temperature is different from (actually is larger than) the equilibrium transition temperature. 
In fact after a quench, the system is almost trapped in those minima and thus displays the aging phenomenon at long 
time. The local minima responsible for the slow dynamics appear at a temperature higher than T^q, thus Td > Teq. 
The activated dynamics which would probably restore the equality Td = Teq is on timescales diverging with N, 




FIG. 7: Spin correlation as a function of r/tw for tw = 15, 20, . . . , 95 (from bottom to top (tmax = WO)), Jb ~ Jh — S = 1 
and T = 0.1. These are the same curves plotted in Fig. ^ but with respect to the variable t jt^. The excellent collapse strongly 
suggests that the function h(t) (present in the dynamical Ansatz) should be equal to t. 



completely unaccessible to our mean- field analysis. 



E. Analysis of the aging regime 



Let us now focus on the aging regime, i.e. t, t' and also t — t' are large but the ratio h{t')/h(t) stays finite. Note 
that within the following asymptotic analysis one cannot find out what is the function h(t) [g[ This is indeed 

an open problem already for classical systems. However, the numerical results, see FIG. 0, suggest that for our model 
h{t) = t. 



Plugging the Ansatz (21a, 21b) into the dynamical equations and after some manipulations similar to p5[ we get: 



Ai?AG(M) 



J:i^{uj = 0)B.AG{f^) + ^^'^{^i)Rst{oj = 0) 



AGf 



: /' ^J:rix)R*^a (-) + f dxlli^{x)KAG (^) + ^W{^ = O)X^g(m) 
Jo A* VA*/ Jo 

+ E^^(/x)ii^^(c. = 0) 



^le^P^) = - X (^IgO^)^Ag(A^) - mAG{n)?RAG(n)) 
^i''{f^)^%KAGM\'KAGit^) 



(28) 

(29) 
(30) 

(31) 



where Kag satisfies the boundary condition Kag{^) 
important to remark that: 



-ig and we have used the notation /i — h{t')/h(t). It is 



1. there is no bath contribution to the aging part of the self-energy. This is natural and it is probably generally true 
since the bath has always its own equilibration timescale, therefore in the aging time-sector (i, t' and t — t' are 
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FIG. 8: TTI part of the spin response function (dots) compared to its expression (continuous line) computed by the QFDR 
from the correlation for Jb = Jh = 1, S = 1, tw > 10 and T = 0.1. 

large but the ratio h[t')/h{t) stays finite) the bath is always already equilibrated and cannot give a non-constant 
aging contribution. 

2. the terms linear in R (and higher) have been neglected in Tij^, whereas the terms quadratic (and higher) have 
been neglected in S^'^. This is due to the fact that they do not give a finite contribution in the aging equations. 

3. As pointed out in the classical case ||] and recently in the quantum case the equations ( p8| , p9| , p0pl| ) are 
re-parameterization invariant. However there is only one function h(t) reached by the system in the long time 
limit. This is a general problem arising in the study of the asymptotic solution of partial and integro-diffcrcntial 
equations, called the matching problem. Until now different techniques arc known and applied to solve this 
problem for partial differential equations but its solution for the dynamical equations arising in the study of 
glassy systems remains an open problem. 

4. The correlation and response functions arc supposed to be related by the GFDR in the aging regime after that the 
exponential e~'^^*~' has been gauged out. In particular the GFDR predicts that Rag{^) — ~ TTTfi-^ ag)' {lA 
(in our notation the GFDR for bosons has a — i/(2Tcfr) instead that 1/Tcff). One can indeed verify that this 
is really a property of the aging equations: (|2^) can be obtained by differentiating (^) and using the GFDR. 
Moreover, we remark that if the bosonic correlation and response functions verify the GFDR so do the spin 
correlation and response functions as in (|2^ ) : Tcff is the effective temperature. 

Evaluating the aging equations in /i = 1 and imposing the existence of an aging solution, i.e. Rag{^) 7^ we obtain 
two matching conditions with the stationary regime pq] : the first one is 

Rst{oo - 0) (32) 

and the second is the same one already obtained from the long-time limit of the stationary equation, i.e. ^ = 0, (eq. 
(p7|)). One can simplify further these two matching equations. Indeed, thanks to the GFDR (p^), one can perform 
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FIG. 9: Spin correlation as a function of t in the fermionic case for — 12, 14, . . . 48 (tmax ~ WO), Jb ~ Jh = 1, S = 0.5 
and T = 0.1. Inset: Spin correlation as a function of t in the bosonic case for the same value of the parameter. This figure 
clearly shows the existence of aging in the bosonic case and its absence in the fermionic one. 



the integrals on the aging functions in A. Hence, using the zero frequency term of (p5|), we get: 



0)5^ 



R 



AG 



Rag 



^k9 



Rsriuj^O) 



T,^R*gj,{u} 



0) 



(33) 
(34) 



Moreover, because Rst is a bosonic response function in a pseudo-equihbrium regime, its zero frequency component is 

real and negative. Therefore is a real number equal to one or minus one that we will note ( in the following 

/I— 1 

and i?sT(w = 0) reads: 



Rsrito = 0) 



(35) 



Plugging this expression onto (Q) and using that — —i-^g'^ we finally obtain the equation for Toff. The aging 
solution corresponds to C = +1 (C = ^1 implies y^/Tcff = 0) and is characterized by the following equation: 



JhQea ( p: 1 



(36) 



where we have replaced — ^qEA- Note that, replacing Toff with T jx (where x is the breakpoint in the one-step 
replica symmetry breaking scheme) this becomes the same equation obtained in |pl| using the marginality condition. 
Indeed (^ 26) with A = and ( ^ , ^ and the boundary condition Kst(j = 0) = —i{2S + 1) + ig are a closed set 
of equations that completely determines Kst, RsR, ^iTcs,qEA- In Appendix [A|, we show that they are completely 
equivalent to the equations studied in |2^, ^ using the marginality prescription within a pure static computation. 
Finally, let us stress that, even if Jb is not present in eq. (pq), T^s depends on Jb because the Edwards- Anderson 
parameter depends on Jb via the eqs. (Ea,E6|) which contains explicitly this coupling constant. 
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IV. ROLE OF THE RELAXATION BATH 



In this section, we briefly discuss the effect of the couphng strength to the relaxation bath Jb- In deriving Eqs 
(p^), we took a generic bath and expand in second order in its coupling constant Jb- Thus the equations we derived 
are a priori only valid in the limit of small Jb- However, we will show that our main equations ( p^ also describe the 
dynamics of a model with finite Jb in a extreme limit. Thus it is legitimate to study them for Jb finite (there is no 
risk of inconsistencies). 

The effect of the bath will of course depend on its precise form. Two types of bath can be considered : baths that 
only couple to the spin, and baths that couple to the boson or the fermions. In this discussion, we will concentrate 
on the first kind, since the spin is the physical object, not the boson. One of the simplest possibility is to couple 
the spin to 2 fermions using the Kondo interaction. In order to take the large-TV limit, we directly introduce the 
SU{N) X SU{N-f) Kondo model, with N-f flavors, defined by g| : 

^ ^ rfjN ^ ' "^i + 51 ^kc\ia'^kia + jr~7= X/ Sia|3cl^pCk' ia (37) 

"^■'V^^ i<j k V ^ l<Q,/3<Ar 

l<i<Af7 KiKNy 
l<a<N 

where c are the bath fermions, efc their kinetic energy and Jb is now the Kondo coupling. In the large- A'^ limit, one 
can still find a closed system of equations, but at the expense of introducing auxiliary f ermionic Gre en functions r 
and k, as explained for example in p3 |. The dynamical equations are similar to (^6|) : (16a 16b, 16d) are the same 



and the bath term in self energies (16i 



16j) are replaced by 



Sfl = ( {k' + li')R* - 2\K\'r] (R*ok - K*r] (38a) 



= -3.(^2K\R\' - K* + K^) ) - (^i?Sr - K*k^ (38b) 



Since the bath has now a proper dynamics, r and k should be computed using new Dyson equations : 



r{t,t')=6{t-t')+ j du<jR{t,u)r{u,t') (39a) 
k{t,t') ^ [ du aR{t,u)k{u,t') + [ du (TK{t,u)r{t' ,u)* (39b) 



and their self-energy reads: 



,Heis 



iJB 
iJB 



RqK + KqR (40a) 



r should also satisfy the boundary condition : 



lim r{t,t') = -i (41) 
t^(t')+ 



It is not difficult to show that in the limit of an infinite number of channel 7 — )■ oo, these equations reduce to (p^). 

However, for finite 7, these equations are much more complex than (|l^), since at low temperature the Kondo 
scale appears and one has to deal with a problem with many different scales. This really increases the difficulty of a 
numerical computation. However, one can extend rather simply the previous analytical study and verify that the same 
dynamical scenario continues to hold. In this paper, we restricted ourselves to study (^6|) as a function of the strength 
of the bath Jb , using Matsubara formalism with marginality condition p2[ . We found that increasing Jb the spin glass 
transition temperature Tg decreases monotonically and as far as we can solve the numerical equations, the transition 
is still second order, given by the condition a; = 1 (x is the value of the breakpoint in the replica formalism). However, 
the decrease of Tg is slow and we could not reach numerically a point where it vanishes. Numerical computation 
can not for the moment decide whether there is a second order phase transition until a quantum critical point at 
finite Jb and T = or the spin glass is not destroyed at zero temperature until Jb = 00. As emphasized in our 
concluding remarks, this situation is disappointing since we would like to study the aging in the vicinity of a non 
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pathological quantum critical point for this model. Our interpretation is that this bath, coupled to the spin directly, 
is not "efhcient" enough and that we probably need to couple to a bath with charge fluctuations by introducing holes 
in the model. Such a doped model is also interesting physically to study the destruction of a quantum spin glass by 
doping. Finally, let us emphasize that solving numerically the model with the Kondo bath ( |37| ) or a more general 
bath may lead to more interesting results, such as an increase of the critical temperature and the Edwards- Anderson 
parameter as the coupling to the bath increases from . 



V. SUMMARY AND DISCUSSION 



In this paper we have studied the out of equilibrium dynamics of the quantum Heisenberg spin glass defined on a 
completely connected lattice and coupled to a spin thermal bath. We have replaced the SU{2) spin symmetry group 
with SU{N) and we have considered the large- iV limit. This has allowed us to have a more tractable model which 
however seems to capture some of the physics of the SU(2) case Thanks to the large- limit we have obtained a 
set of closed integro-differential equations on the correlation and response functions. By the analytical study and the 
numerical integration of these equations we have fully analyzed the real time (dissipative) dynamics of the mean field 
quantum Heisenberg spin glass model in the large- A^ limit. We have considered a particular type of initial condition 
which corresponds to the physical situation in which, at i = 0~ the system is at equilibrium at infinite temperature 
and at i = 0"^ becomes coupled to thermal bath in equilibrium at temperature T. This corresponds to an extremely 
fast quench from very high temperature. Depending on the value of T, the system has a very different long-time 
behavior. 

At high temperature the system relaxes, after a finite equilibration time, inside the paramagnetic state. In this 
stationary regime the system is at equilibrium and the fiuctuation dissipation relation holds. When the system is 
quenched below a certain critical temperature Td, which depends on the values of the spin and the system-bath 
coupling, it never reaches an equilibrium regime. At large times two time-sectors can be identified for the behavior 
of the correlation Cs{t,t') and response Rs{t,t'). When t and t' are very large, but their difference remains of the 
order of one, the systems reaches a pseudo-equilibrium regime in which the QFDR is verified. However on a larger 
timescale, diverging with the age of the system (t — i' cx t,t'), there is a secondary relaxation called aging. In this 
regime the quantum fluctuation dissipation relation is violated and the correlation and the response are related by a 
generalization of the classical fluctuation dissipation relation characterized by an effective temperature different from 
the bath temperature. 

Moreover we have also studied the role of the bath. First, we have taken a linear coupling of the spin to the bath, 
we have developed to the second order in the coupling constant and integrated out the bath spins. In this way we 
have found a generalization of the Feynman- Vernon influence functional p5[ | for spins in which the properties of the 
bath enters only through its susceptibility. All the numerical study has been done in this case. However, we have 
also considered a more general type of bath and we have shown that a "simple" one turns out to be a Kondo bath. 
We have extended the analytical study to this case and shown that the previous dynamical scenario continues to 
hold. Furthermore we have unveiled that the way we have followed previously to treat the system-bath coupling can 
be recover as a limiting case of a Kondo Bath. Finally, we have also verified numerically that, as far as we can go 
increasing the coupling to the bath in (^, the dynamical transition remains of second order (by this we means that 
the asymptotic dynamical energy is continuous) and the critical temperature does not vanish. 

The most striking features of the low temperature out of equilibrium dynamics are the aging phenomenon and the 
generalization of the fluctuation dissipation relation out of equilibrium. It has been shown for spherical spins p5| ] and 
for rotors |2^, Q that a generalization of the classical fluctuation dissipation relation holds in the aging regime. In 
this paper we have shown that this is the case also for models with a non trivial spin algebra. These results seem 
to suggest that, except for the renormalization of the coefficients of the dynamical equations, the aging regime is not 
affected by quantum fluctuations and the aging systems behaves classically in their slow evolution. But is this always 
true? Is it not possible to find "a quantum system which ages coherently" ? Since in general, the decoherence time is 
finite and the aging regime takes place in the large time limit, a classical aging regime is always expected to set in at 
large enough time. However there is an important case in which this naive argument may fail. Near a quantum critical 
point the decoherence time diverges, therefore it could be possible that at very large times (larger than the time on 
which the system enters in the asymptotic regime and than the characteristic timescale of the TTI- regime), but still 
lower than the decoherence time, the system ages coherently. We could not address this very interesting question for 
the quantum Heisenberg spin glass analyzed in this paper : the technical reason is that its quantum critical point is 
rather pathological since it corresponds to a vanishing spin size. Hence, another type of model with a less singular 
quantum critical point has to be studied. Work is in progress in this direction p6|. 
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APPENDIX A: FROM REAL TIME TO IMAGINARY TIME 

In this appendix, we give explicit formulas for doing the Wick rotation to imaginary time in equilibrium. Let us 
define : 

^^J^*^- lG__(t)fort<0 



We note that this function has a simple expression in terms of the spectral density (using equilibrium FDT) : 



i J dep{e)hB{£)e for bosons 
—i J dep{e)hp{e)e''^'^* for fermions 



[G]{t) = 



—i j dep{e)nB{f)e.'^^^ for bosons 
i J dep{e)np{e)e'^'^* for fermions 



where ^^(e) and np{e) are the Fermi and Bose function respectively and we use the notation f{x) = f{~x). Thus 
[G] is analytic in t, and we have the relation : 

[G]{~iT) = iG{t) 0<t<(3 (A2) 

\G](-iT) = / ^^^^ 0<T < (3 for bosons 

^ I -iG{P~T) {)<T<(3 for fermions ^ ' 
where the Matsubara Green function is defined by |47j : 

G{t) = - {Tb{T)h\Q)) ^ j de p{e)hB{e)e~" for bosons (A4) 

G(t) EE - {Tf{T)f{0)) = - /" de p{€)nF{e)e-''' for fermions (A5) 



The same formula also holds for the self-energy. 

Using this result, we find the imaginary time equations in the paramagnetic state : 

{G-^){iiyn) = w„ + A - I](ii/„) (A6a) 

E(r) = G(r) (^J%G{r)Gi^r) + J|x"(r)) (A6b) 

G(t = 0-) = -S' (A6c) 

They are a slight generalization of Eq. (5) of |^ , including the bath. 

Similarly, in the glassy phase, using the same technique, we find Eqs. (31) of |Q with 8 = which corresponds 
to the marginality criterion, as explained in . 

APPENDIX B: NUMERICAL SOLUTION 

In this appendix, we provide some details about the numerical solution of our main equations (|l^). In order to 
compute R and K on the domain O < i' < i, we use the causality of the equations (|l^) : in order to compute the 
function in (t, t') we only need the knowledge of the functions at previous times, so we can construct the functions 
step by step in time along the t direction. This structure of the equations is general for classical or quantum spin 
glass dynamical problems. (See e.g. |^). We have to solve a set of coupled differential equations in t. However, in 
this problem the situation is more complicated, since we have first to compute an unphysical bosonic function, which 
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oscillates a lot. A naive algorithm is to compute the derivatives a each point {t,t') for a fixed t, and extrapolate 
using a first order Taylor expansion. However for numerical integration of ordinary differential equations (ODE), this 
method is not recommended (See e.g. [|8|), since one needs a very tiny mesh size to obtain accurate result. In our 
case, we found that this simple algorithm does not give any good result for a reasonable computational cost, contrary 
to simpler models studied previously (e.g. classical p— spins models). 

Hence, we used a modified procedure, inspired by the Stoer-Burlish algorithm for (ODE) : let us assume that we 
have computed the functions until time t and we want them at time t -\- 6 where S is our mesh size. We cut this 
step into N parts, and compute the functions for t + iS/N for all t' and 1 < i < N ^ using the modified midpoint 
method |Q. We then obtain the functions at t + S, for various N and all t' and we extrapolate the result to — > oo. 
Typically, we use 3 or 4 values of N among {4, 8, 16, 32}. The integrals are computed using either a trapezoidal or a 
Simpson formula. It is important to notice that the structure of the equations (|l^) implies that we do not need to 
keep the intermediate point after the t + S have been computed. As explained in the text, the dynamical equations 
conserve the constraint, which is automatically satisfied in the time evolution : it is also important for the stability 
of the algorithm that its discrete implementation of Eqs. (16) respects this conservation exactly. 
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